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Abstract: In this brief report, we analyzed ecological correlates of risk of extinction for mammals endemic to China using 
phylogenetic eigenvector methods to control for the effect of phylogenetic inertia. Extinction risks were based on the International 
Union for Conservation of Nature (IUCN) Red List and ecological explanatory attributes that include range size and climatic 
variables. When the effect of phylogenetic inertia were controlled, climate became the best predictor for quantifying and evaluating 
extinction risks of endemic mammals in China, accounting for 13% of the total variation. Range size seems to play a trivial role, 
explaining ~1% of total variation; however, when non-phylogenetic variation partitioning analysis was done, the role of range size 
then explained 7.4% of total variation. Consequently, phylogenetic inertia plays a substantial role in increasing the explanatory power 
of range size on the extinction risks of mammals endemic to China. Limitations of the present study are discussed, with a focus on 


under-represented sampling of endemic mammalian species. 
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Species’ extinction risk is driven by multiple ecolo- 
gical factors (Cardillo et al, 2005, 2008), including small 
population size (Legendre et al, 2008), small range size 
(Harris & Pimm, 2008), large body size (Cardillo & 
Bromham, 2001) and degrading habitat conditions 
(Halley & Iwasa, 2011). One key focus of current 
macroecological studies is understanding the ecological 
correlates of species’ extinction risks, in order to both 
better understand change in demographics and to 
implement more effective conservation measures 
(Cooper et al, 2008; Keane et al, 2005; Reynolds et al, 
2005). 

Global patterns of the extinction risks of mammals 
have been well quantified (Cardillo et al, 2004, 2005, 
2006; Jones et al, 2009), but it is unknown whether 
the relevant ecological determinants attributed to 
extinction risks of mammals on a global scale can be 
applied at regional or local scales. The central goal of 
the present study is to evaluate drivers of extinction 
risks of mammals at the regional scale, specifically 
examining the ecological causes of extinction risk for 
mammals endemic to China and accounting for 
phylogenetic inertia (Carrascal et al, 2008). There are 
over 200 known mammal species endemic to China, 
with the full distributional ranges only limited to the 
terrestrial boundaries of modern China. However, 
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because of limited data access, only a subset of 
endemic mammals has been included in the present 
analysis, though I deal with the corresponding pote- 
ntial constraints. 


MATERIALS AND METHODS 


Data sets 

Distributional records of endemic mammals were 
derived from the China Species Information Service 
(http://www.baohu.org/) and literature related to 
mammalian fauna of China (Smith & Xie, 2008; Wang, 
2003). For the selected species, phylogeny was recons- 
tructed from the previously established meta-phylogeny 
of global mammals (Bininda-Emonds et al, 2007). The 
tree for the subsequent analyses is presented in Figure 1. 
The threatened status of each species was obtained from 
the International Union for Conservation of Nature 
(IUCN; http://www.iucn.org/). 
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Eozapus setchuanus 
Apodemus semotus 
Apodemus chevrieri 
Niviventer andersoni 
Niviventer culturatus 
Cansumys canus 
Eothenomys chinensis 
Eothenomys custos 
Eothenomys olitor 
Eothenomys proditor 
Lasiopodomys fuscus 
| Proedromys bedfordi 
Volemys millicens 
Volemys musseri 
Brachiones przewalskii 
Meriones chengi 
Sciurotamias davidianus 
Sciurotamias forresti 
Petaurista alborufus 
Petaurista xanthotis 
Trogopterus xanthipes 
Aeretes melanopterus 
Lepus comus 

Lepus hainanus 
Ochotona erythrotis 
Ochotona muliensis 
Ochotona himalayana 
Ochotona iliensis 
Ochotona cansus 
Ochotona gaoligongensis 
Ochotona thomasi 
Macaca thibetana 
Macaca cyclopis 
Pseudois schaeferi 
Procapra przewalskii 
Elaphurus davidianus 
Muntiacus crinifrons 
Lipotes vexillifer 
Ailuropoda melanoleuca 
Rhinolophus osgoodi 
Rhinolophus rex 
Plecotus taivanus 
Murina puta 

Uropsilus andersoni 
Uropsilus investigator 
Uropsilus soricipes 
Mogera insularis 
Euroscaptor longirostris 
Scaptochirus moschatus 
Scapanulus oweni 
Mesechinus hughi 
Sorex excelsus 
Blarinella quadraticauda 
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Figure 1 The phylogeny of 53 mammals endemic to China 
(redrawn from Binind-Emonds et al, 2007) 


Ecological variables 

Range size for each endemic mammal was calculated 
using digital distribution range maps (http://www. Iuc- 
nredlist.org/technical-documents/spatial-data). Climatic 
variables included precipitation, minimal temperature, 
maximal temperature, mean temperature, evaporation, 
humidity and solar radiation. These data are interpol- 
ations of observed data collected from 1950-2000. Data 
were calculated and exacted from grid cells with the 
presence of the species using the WorldClim database 
(http://www.worldclim.org). Data are available from the 
author upon request. 


Data analyses 
I used phylogenetic eigenvector regression (PVR) 
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(Carrascal et al, 2008; Diniz-Filho et al, 1998, 2011; 
Kuhn et al, 2009; Morales-Castilla et al, 2012; Seger 
et al, 2013) to quantify and account for phylogenetic 
signals inherited in tip species due to non-independent 
history caused by cladogenesis and anagenesis events 
(Inglis, 1988). The core of the PVR method is to 
construct a matrix of pairwise phylogenetic distances 
among concerned species and a principal coordinates 
analysis (PCoA) is performed on this matrix. The 
matrix of phylogenetic distances for the endemic 
mammals of China is available from the author upon 
request. The most important axes accounting for most of 
total variance were retained to represent the major 
components of phylogenetic signal for subsequent 
analyses (Carrascal et al, 2008; Diniz-Filho, et al, 2012c). 

To remove phylogenetic autocorrelation within the 
climatic variables and life-history traits for endemic 
mammals, all the explanatory variables were regressed 
onto the selected eigenvectors so as to obtain residuals to 
study the correlations between explanatory variables and 
extinction risk. 

Variation partitioning followed the methods of 
previous studies (Legendre & Legendre, 1998; 
Carrascal et al, 2008). In detail, climatic variables 
were used as group 4, while life-history traits formed 
group B. During the multiple regression analysis when 
extinction risks of species served as the response 
variable, the difference between the variation V4, 
explained by A+B and the variation V4 explained by A 
became the independent contribution of group B for 
explaining the total variation of extinction risk. 
Similarly, the independent contribution of group A for 
explaining the total variation of extinction risks of 
species was given by V4,—Vz. Unexplained variation 
inside the response variable was 1—V4;. 

As a comparison, variation partitioning without cont- 
rolling phylogenetic inertia was also performed to 
evaluate the relative influence of phylogenetic inertia. 
The variation partitioning procedure was identical to the 
above-mentioned method, except that all variables were 
not phylogenetically corrected. 


RESULTS 


When phylogenetic inertia was controlled using 
PVR, there was still a large fraction of unexplained 
variation for the extinction risks of endemic mammals of 
China, with only 13% of total variation being explained 
by the ecological variables included here (Figure 2). The 
influential role of climate was identified, which explains 
approximately 13% of the variance independently. As 
such, almost all the fraction of explained variation was 
attributed to climate. By contrast, the role of range size 
was not important, because only 1% of total variation 
was explained by this trait (Figure 2). 
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Figure 2 Partitioning of variation attributed to different 
explanatory variable groups after controlling for 
phylogenetic inertia effects 


When variation partitioning was performed without 
controlling for phylogenetic inertia (Figure 3), the results 
were similar to those resulting from phylogenetic 
adjustment (Figure 2), but the role of range size in 
extinction risk became more important, accounting for 
7.4% of the total variation in extinction risk and climate 
explained 11.2% of total variation (Figure 3). 
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Figure 3 Partitioning of variation attributed to different 
explanatory variable groups without controlling for 
phylogenetic inertia effects 


As a consequence, the influence of phylogenetic 
inertia could be determined via comparison of these two 
methods (contrast Figure 2 and Figure 3). The variation 
partitioning results for the situation with phylogenetic 
adjustment would partially eliminate the influence of 


range size when explaining the extinction risks of species. 


DISCUSSION 


There are a suite of statistical methods to detect and 
remove phylogenetic inertia effects (Blomberg et al, 2003, 
2012; Diniz-Filho et al, 1998, 2012b; Pagel, 1999), but 
many are strongly correlated and robust for handling 
alternative evolutionary models (Diniz-Filho et al, 2012b; 
Seger et al, 2013). The merit of the PVR method and its 
extension (Diniz-Filho et al, 2012a) is that it is able to cap- 
ture nonlinear components of phylogenetic signals and 
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elegantly evaluate the deviation of Brownian motion for the 
evolution of those examined traits (Diniz-Filho et al, 2012a). 

Previous studies found that range size is an important 
attribute that can predict extinction risk over different 
taxonomic groups (including global mammals) (Cardillo et 
al, 2005; Cooper et al, 2008; Harris & Pimm, 2008; Hanna 
& Cardillo, 2013; He, 2012). However, for the endemic 
mammals of China, threatened status is principally 
predicted by climate, and not life-history or range size (Figs. 
2 and 3). One reason for the weak explanatory power of 
range size in my findings is that the distribution information 
for endemic taxa in China may be incomprehensive, 
reducing the explanatory power of range size on structuring 
extinction risk. Furthermore, limited phylogenetic sampling 
of endemic species may have contribute to the reduced 
explanatory power of range size. 

Despite interesting findings, there are a few 
shortcomings to note. First, the weak phylogenetic inertia 
signal detected may be due to the fact that many variables 
influence species’ extinction risk besides or in addition to 
climate and distribution; for example, human disturbance 
(Kong et al, 2013) and environmental pollution. A lack of 
comprehensive and comparable data on such topics for use 
in examining extinction risk makes for a less convincing 
or robust analysis. Moreover, since I did not include non- 
endemic species in the analysis, the phylogenetic signal is 
weakened because the role of evolution is weak when 
influencing species’ range sizes, due to limited 
evolutionary divergence amongst endemic taxa. 

A further potential limitation of the present study is 
the comparatively small sample of 53 endemic mammals, 
conspicuously less than the total number of endemic 
mammals in China. As a consequence, these findings 
may not be robust enough to represent the true ecological 
correlates of extinction risk for this group. The power of 
phylogenetic signals is also reduced due to this limited 
sampling. Unfortunately, there is no clear remedy to this 
shortcoming, as only 53 endemic mammalian species 
were able to be included due to data, e.g., global 
mammalian phylogeny does not contain many mammals 
endemic to China (Bininda-Emonds et al, 2007), some of 
which may be yet identified as new endemic species or 
become a synonym/subspecies of a wide-ranging species 
(Li et al, 2006). Further research in the area may allow for a 
more detailed replication of this analysis at a later date. 
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